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ABSTRACT 


Intermit  tency  measurements  have  been  performed  within  the  first  one 
hundred  virtual  diameters  of  an  axi-symmetr ic  compressible  wake.  The 
primary  measurements,  performed  with  the  hot-wire  anemometer,  concerned 
the  intermittency  factor  and  the  frequency  of  zero  occurrences.  The 
intermittent  flow  covered  most  of  the  wake  profile  but  for  a  narrow 
region  about  the  axis,  and  was  found  to  be  normally  distributed  around 
the  average  front  position.  The  latter,  along  with  the  extent  of  the 
front  standard  deviation,  was  found  to  agree  numerically  with  expecta¬ 
tions  based  on  Low  speed  wakes  and  to  grow  as  the  1/3-power  of  axial 
distance.  A  weak  periodicity  of  the  front  was  detected  at  a  wavelength 
about  nine  times  greater  than  the  longitudinal  scale  of  turbulent  velocity 
fluctuations.  This  periodicity  affects  the  turbulent,  spectra. 
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FOREWORD 

This  report  is  the  third  in  a  series  describing  the  properties  of  free 
compres s ib Le  turbuLent  flows,  including  turbulent  plasmas,  stemming  from 
a  continuing  research  program  on  hypersonic  wakes  at  Philco-Ford 
Corporation.  The  present  work  describes  the  intermittency  character¬ 
istics  of  an  axi-symme tr ic  compressible  wake  whose  mean  and  turbulence 
properties  have  been  described  in  the  two  reports  preceding  it,  Mean 
Flow  Measurements  in  an  Axisymmetric  Compressible  Turbulent  Wake" 

(No.  U-3978,  1  March  1967)  and  "Turbulence  Measurements  in  an  Axi~ 
Symmetric  Compressible  Wake"  (No.  UG-4118,  1  August  1967).  The 
corresponding  properties  of  a  two-dimensional  wake  will  be  reported  soon. 

This  research  has  been  supported  by  the  Advanced  Research  Projects  Agency 
and  monitored  by  the  Space  and  Missile  Systems  Organization,  USAF .  The 
author  acknowledges  with  gratitude  the  continuing  encouragement  of  the 
Philco-Ford  Corporation  and  the  expert  assistance  of  Mr.  Lee  Von  Seggern 
with  the  measurements  and  the  instrumentation. 
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NOMENCLATURE 

coefficient  of  polynomial  curve-fitting  the  intermittency 
factor 

coefficient  of  polynomial  curve-fitting  the  zero  frequency 

(square  of)  virtual  wake  diameter 

frequency;  also,  random  function  of  Appendix  B 

frequency  at  front  wavelength 

jump  across  the  front  (Appendix  B) 

transverse  scale  of  wake 

wake  half-radius 

local  Mach  number 

frequency  of  zero  occurrences 

value  of  N  at  Y 

maximum  of  N  at  each  axial  position 
probability 

turbulent  Reynolds  number 
Local  temperature 
non-dimensional  temperature 
Local  axial  velocity 
non-dimensional  axial  velocity 
free  stream  axial  velocity 
axial  velocity  at  Y 

computer  program  for  data  reduction  (Appendix  A) 
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Axial  coordinate 

non-dimensional  axial  coordinate 
axial  position  at  which  data  were  taken 
radial  coordinate 
radial  position  of  the  front 

radial  position  at  which  radial  derivative  of  y  is  zero 

radial  position  at  which  N  is  maximum 

intermittency  factor;  also,  ratio  of  specific  heats 

r.m,s.  fluctuation 

non-dimensional  radial  coordinate 

front  autocorrelation  microscale 

average  wavelength  of  front 

longitudinal  macroscale  of  velocity  fluctuations 
longitudinal  macroscale  of  temperature  fluctuations 
local  density 

standard  deviation  of  front 
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1 .  INTRODUCTION 

It  is  well  known  chat  turbulent  flows  spreading  into  ambient,  quiescent 
fluid  are  separated  from  the  latter  by  a  sharp,  irregular  front  or 
interface.  This  front  changes  shape  continuously  with  time,  but  on  the 
average  it  advances  into  the  ambient  fluid  engulfing  the  latter  and 
making  it  turbulent.  In  shear  turbulence  this  motion  is  accompanied  by 
an  overall  movement  of  the  front  along  the  direction  of  the  shearing 
motion  as  well.  Thus,  there  are  at  least  three  practical  objectives 
motivating  the  study  of  this  front.  First,  its  statistical  properties 
are  necessary  for  the  quasi-steady  geometrical  description  of  shear 
zones  such  as  boundary  layers,  jets  and  wakes.  Secondly,  the  correct 
measurement  of  the  net  turbulence  properties  within  the  front  requires 
the  recognition  and  separation  of  those  effects  on  the  turbulence  signal 
which  are  due  to  the  intermittent  nature  of  the  front  itself.  Finally, 
the  more  detailed  the  understanding  of  ambient  fluid  enguifment  and  con¬ 
version  is,  the  better  will  be  the  prediction  of  the  bulk  properties  of 
the  shear  flow,  such  as  thickness  and  lateral  expansion.  These  objectives 
are  particularly  important  for  the  hypersonic  wake  vhere  the  shape  of  the 
front  envelope,  its  motion  and  its  effect  on  the  turbulence  within  it 

enter  directly  and  quantitatively  in  computations  of  radar  backscatter. 

« 

This  report  deals  primarily  with  the  geometrical  features  of  the  front 
and  more  specifically  with  the  most  probable  front  position,  the  standard 
deviation  (a  measure  of  the  "wrinkle"  amplitude),  and  the  front  micro¬ 
scale  (a  measure  of  the  front  wavelength).  A  considerable  investigation 
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found  can  only  be  related  to  gross  flight  parameters  (e.g.,  body  shape 
and  size)  rather  than  the  more  meaningful  integral  wake  properties  (e.g., 
the  drag,  transverse  scale,  etc.)  obtainable  by  a  complete  mapping  of  the 
flow  field  such  as  is  possible  in  wind  tunnels.  The  latter  prerequisite 
has  already  been  fulfilled  in  the  present  case,  since  complete  flow  maps 
of  the  axi-symmetric  wake  at  hand  were  already  available  (References  9 
and  10).  Using  this  well-probed  wake  flow,  two  additional  measurements 
were  made:  the  intermittency  factor  was  recorded  with  a  special  elec¬ 
tronic  circuit,  and  the  frequency  of  "zero-occurrences"  was  also  measured. 
Sufficient  information  was  thus  gathered  to  deduce  (1)  the  position  of 
the  turbulent  front,  (2)  the  front  standard  deviation,  and  (3)  the  front 
microscale.  Further  use  of  these  data  can  and  will  be  made  in  deducing 
further  turbulence  properties  of  the  wake. 

2.  WAKE  MODEL  AND  FLOW  FIELD 

The  wind  tunnel,  asixymmetric  model  and  wake  flow  characteristics  have 
already  been  described  in  References  9  and  10.  The  reader  should  consult 
these  references  for  details  and  especially  for  th*at  nomenclature  not 
peculiar  to  the  intermittency  measurements  themselves.  A  photographic 
and  schematic  view  of  the  model  appears  on  Figure  1. 

3.  PROCEDURE  AND  SPECIAL  INSTRUMENTATION 

The  basic  tool  for  intermittency  measurements  is  a  method  capable  of 
recording  the  time  the  hot-wire  is  immersed  in  the  turbulent  flow  versus 
the  time  it  is  not.  The  simplest  way  to  do  this  is  to  obtain  a  continuous 
direct  record  of  the  interface  (e.g.,  hot-wire  output  versus  time  or  a 
series  of  wake  photographs)  and  to  subject  this  record  to  statistical 
analysis;  this  method  is  in  general  awkward  and  inaccurate.  A  better  way 
is  to  perform  the  statistics  in  situ  by  means  of  an  electronic  circuit 
such  as  used  by  Corrsin  and  Kistler  (Reference  1)  and  Bradbury 
(Reference  11).  The  device  described  below  is  basically  their  design 
adapted  to  higher  frequencies. 
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FIGURE  1.  TOP:  PHOTOGRAPH  OF  TEST  SECTION  SHOWING  AXI SYMMETRIC  MODEL 
AT  LEFT,  PROBE  AT  LEFT  AND  PROBE  ACTUATOR  (TOP  LEFT). 

BOTTOM :  DETAILS  OF  THE  AXISYMMETRIC  MODEL. 


of  these  features  has  been  performed  it  low  speeds  with  such  flows  as  jets, 
wakes,  and  boundary  layers:  compendia  of  these  studies  have  been  presented 
by  Corrsin  and  Kistler  (Reference  l),  and  more  recently  by  Townsend 
(Reference  2).  These  compendia  show  that  the  qualitative  features  of 
the  front  irregularity  follow  some  general  rules  irrespective  of  flow 
geometry;  on  the  other  hand,  quantitative  differences,  such  as  the  ratio 
of  the  standard  deviation  a  to  the  average  front  position  Y,  have  beer 
noted  among  boundary  layers,  jets,  and  wakes.  Furthermore,  corresponding 
front  characteristics  for  the  axisymmetric  wake  have  not  been  reported  to 
date  in  comparable  thoroughness  and  detail,  although  some  tentative 
intermi ttency  observations  in  sphere  wakes  have  been  reported  by 
Baldwin  (Reference  3)  and  Gibson  (Reference  4).  A  curious  phenomenon 
recorded  in  the  course  of  these  studies  is  the  appearance  of  rather 
extensive  stretches  of  irrotational  flow  on  the  wake  axis  where  fully 
turbulent  flow  is  usually  noted  in  the  results  of  other  experiments  of 
this  type  (see  Reference  l) .  Attempts  to  examine  this  phenomenon  in  the 
supersonic  wake  were  made  in  the  present  experiment. 

The  main  objective  of  this  portion  of  the  present  work,  however,  was  to 
study  the  effect  of  compressibility.  High-speed  measurements  of  the  front 
structure  are  virtually  nonexistent,  with  the  exception  cf  the  statistical 
interpretation  of  ballistic-range  photographs  of  wakes  carried  out  to 
some  detail  by  Schapker,  Levensteins  and  Krumins,  and  Knystautas 
(References  5  through  7)  respectively.  This  optical  method  has  revealed 
certain  important  similarities  between  low-  and  high-speed  front  behavior 
as  regards,  for  example,  the  distribution  of  intermi ttency  factor  and 
the  front  microscale.*  Unfortunately  the  intermittent  properties  so 

There  are  several  serious  objections  to  the  measurement  o^  intermittency 
by  the  processing  of  photographs.  For  example,  proper  statistics  of  a 
stationary  ensemble  can  be  performed  only  if  there  are  no  gradients  of  the 

mean  or  root-mean-square  values  in  the  ensemble- i .e . ,  if,  for  instance,  the 
wake  thickness  is  constant  with  distance;  this  is  not  so  restrictive  for 

small-scale  turbulence  studies3  but  for  the  interface,  where  the  scales  are 
large,  it  means  that  proper  statistics  can  be  performed  from  photographs  on 
past  some  thousands  of  diameters.  There  are  also  the  usual  observational 
difficulties:  prominences  and  ridges  can  be  seen  from  the  photographs,  but 
''valleys"  cannot. 


An  intermittoncy  meter  for  measuring  the  intermi ttency  factor  of  turbulent 
wakes  directly  from  a  hot-wire  anemometer  signal  was  fabricated  for  use 
in  studying  turbulent  wake  structure.  The  instrument  utilizes  seven 
Burr-Brown  operational  amplifiers  and  associated  feed-back  networks  for 
signal  conditioning.  A  signal  level  of  up  to  10  volts  peak-to-peak  can 
be  handled  effectively  by  this  system  making  it  compatible  with  the 
Transmetrics  No.  6401-a  hot-wire  anemometer  amplifier.  All  Burr-Brown 
Amplifiers  used  in  the  system  to  condition  the  ac  signal  are  100  volts/ 
usee  slewing  rate  providing  an  extremely  fast  response  instrument  for 
measuring  the  intermit tency .  This  device,  illustrated  in  Figures  2  and  3, 
also  included  a  0-to-l-volt  dc  output  directly  proportional  to  the 
intermit tency  factor,  which  could  thus  be  read  directly  on  the  face  of  a 
panel  meter  or  used  to  drive  a  dc  recorder.  It  was  initially  hoped  that 
the  wire  could  be  traversed  continuously  through  the  wake,  and  that  the 
dc  output  noted  directly  on  the  recorder.  As  other  investigators 
(Reference  1)  also  found,  however,  the  trigger  level  and  other  circuit 
controls  required  adjustment  from  point  to  point,  so  that  the  intermi t tency 
measurement  was  not  quite  that  simple.  Nevertheless,  the  circuit  performed 
satisfactorily  despite  the  high  frequency  response  required. 

For  these  measurements,  an  oven-calibrated  0.00005-inch  wire  is  used.  This 
wire  was  processed  as  indicated  in  Reference  9  but  was  not  flow-calibrated 
since  the  on-off  process  of  intermit tency  requires  no  quantitative  knowledge 
of  the  wire  properties  as  such.  The  wire  current  was  first  adjusted  as 
the  maximum  safe  value  so  that  no  melting  would  occur  anywhere  across  the 
wake  and  its  time  constant  adjusted  to  the  proper  value;  it  was  then 
lowered  at  a  point  in  the  wake  and,  after  suitable  adjustment  of  the  circuit 
cor  :rols,  the  intermittency  factor  (or  the  frequency  of  zeros)  was  noted. 

The  block  diagram  employed  is  shown  in  Figure  4. 
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FIGURE  2.  PHOTOGRAPH  OF  INTERMITTENT  CIRCUIT  (CENTER) 
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FIGURE  3.  ELECTRONIC  CIRCUIT  OF  THE  INTERMITTENCY  FACTOR 
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FIGURE  4.  DIAGRAM  OF  ELECTRONICS  USED  FOR 
INTERMITTENCY  MEASUREMENT 
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3.1  PROCEDURE  FOR  INTERMITTENCY  MEAS UREMEl.'T 

The  intermittent  factor  0  <  y  <  1  was  measured  at  each  of  about  22 
radial  position  and  at  each  of  the  14  axial  positions  (XSTATIONS  0 
through  13)  at  which  all  measurements  with  the  axisymmetric  wake  were 
made.  At  each  point  the  signal  was  first  observed  on  the  Tektronix  type 
549  storage  oscilloscope  for  the  purpose  of  making  the  required  adjust¬ 
ments  in  signal  amplification  and  especially  in  the  trigger  level.  The 
trace  storage  technique,  a  new  development  in  osc i  1  lography ,  proved  an 
invaljable  aid  in  recording  y  correccly.  R>  freezing  the  trace  on  the 
scope  screen  for  an  arbitrary  controllable  period  (from  a  fraction  of  a 
second  to  one  hour)  the  operator  can  determine  much  better  than  hitherto 
possible  whether  the  Schmidt  circuit  triggers  properly.  Typical  dual 
traces  thus  obtained  are  shown  in  Figure  5.  In  these,  the  up  position 
of  the  trigger  output  trace  marks  the  time  the  wire  is  within  the 
turbulent  fluid;  the  ratio  of  this  time  to  the  down  time  is  given  by 
the  circuit  directly  as  the  fraction  of  one  volt  dc  (e.g.,  0.68  volt 
denotes  V  ■  0 .68)  . 

3.2  PROCEDURE  FOR  MEASURING  THE  ZERO  OCCURRENCES 

The  rectangular-pulse  shape  of  the  trigger  output  shown  on  Figure  5 
affords  a  relatively  simple  way  by  which  the  frequency  of  the  wire 
crossings  across  the  wake  front  can  be  measured.  This  was  done  by 
directing  the  trigger  output  signal  into  a  counter  whose  own  trigger 
level  was  set  below  l  volt  dc .  The  crossing  frequency  (also  called  the 
frequency  of  zero  occurrences  or  zeros)  was  then  r ead  directly  by  the 
counter.  This  measurement  was  also  made  at  each  of  the  14  XSTATIONS, 
although  at  fewer  radial  positions  for  each.  Reading  the  zero  frequencies 
required,  of  course,  the  same  adjustments  of  the  intermittency  circuit 
and  amplifier  as  during  the  measurement  of  y. 
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4.  DEFINITIONS  OF  FRONT  PROPERTIES 

A  thorough  discussion  of  the  processing  and  interpretation  of 
intermi t tency  measurements,  together  with  the  necossar>  background 
from  the  theory  of  stationary  random  variables,  has  been  given  by 
Corrsin  and  Kistler  (Reference  1).  The  formulas  necessary  to  process 
the  data  as  obtained  above  will  be  given  below. 


Because  the  probability  P  of  the  existence  of  the  turbulent  front  is 

P  »  dy/dY  (l) 


where  Y  is  the  radial  coordinate,  it  follows  that  the  most  probable 
position  Y  of  the  front  is 


00 


(2) 


The  quantity  Y  is  very  crucial  to  the  geometric  description  of  a 
turbulent  wake  and  is,  in  tact,  the  only  rational  definition  of  the 
surface  bounding  the  turbulent  fluid  in  the  average  or  stead\-state 

sense . 


Another  important  definition  is  the  mean  standard  deviation  of  the  front 
corrugations 


a  =  [f  (Y-Y)2  ^  d  (Y-Y)] 


(3) 


2  2 

where  the  lower  limit  of  zero  is  proper  only  if  dv / dY  and  dy/dY  are 

also  zero  at  Y  *  0.  In  addition  to  Y  and  cr,  a  third  quantity  of 

geometrical  importance  is  the  autocorrelation  microscale  of  the  front 

U 


<*> 
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where  U__  is  the  mean  flow  velocity  at  the  front  location  (that  is, 

Y 

U__  =  U(Y))  and  whe’-e  N  is  the  average  number  of  Eeros  at  Y  a  Y .  Again 
tKe  latter  definition  is  appropriate  to  a  normal  (Gaussian)  distribution 
of  y  and  its  moments. 

5.  RESULTS 

5.1  INTERMITTENCY  MEASUREMENT 

Figure  6  shows  the  distribution  with  radius  of  the  intermittency  factor  y 
at  each  of  the  14  axial  positions  examined.  Nearer  the  body  the  region 
of  fully  turbulent  fluid  (y  =  1)  is  noticeably  broader  around  the  axis, 
but  further  downstream  \  departs  from  unity  almost  immediately  off  the 
axis.  It  appears,  therefore,  that  in  the  se 1 f -preserving  region 
(established  in  Reference  9  to  lie  beyond  x  =  40)  the  fully  turbulent 
fluid  is  confined  to  the  very  immediate  vicinity  of  the  wake  axis;  this 
is  in  general  agreement  with  the  two-dimensional  wake  at  low  speeds 
(Reference  12).  The  reader  should  not  misinterpret,  of  course,  the  fully 
turbulent  region  as  enclosed  by  a  smooth  surface  of  vanishingly  small 
radius;  statistically,  this  region  can  at  times  be  of  diameter  comparable 
to  the  wake  diameter.  All  Figure  6  implies  is  that  the  occurrence,  if  any, 
of  laminar  spots  on  the  axis  is  highly  improbable,  in  noteworthy  contrast 
to  the  extent  to  which  such  spots  appeared  in  the  sphere  wakes  of  Gibson 
(Reference  3)  and  Baldwin  (Reference  4).  In  fact,  a  separate  experiment 
was  performed  wherein  154  oscillograms  of  the  turbulence  on  the  axis  were 
obtained,  of  which  n~"  one  showed  evidence  of  laminar  (external)  flow. 

Each  of  these  was  of  200  microseconds  duration  which,  at  a  wake  thickness 
of  order  1  centimeter  and  a  stream  speed  of  about  60,000  centimeters  per 
second,  confirms  conclusively  the  high  improbability  of  y  <  1  on  the  axis. 

Figu  re  7  shows  (in  more  detail  than  Figure  6)  a  comparison  (at  a  typical 
x)  of  the  distribution  in  y  and  in  the  velocity  u  and  temperature  T,  the 
latter  two  as  obtained  in  Reference  9.  It  is  significant  that  turbulent 
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RADIAL  DISTRIBUTION  OF  INTERMITTENCY  FACTOR, 
WITH  VELOCITY  AND  TEMPERATURE  PROFILES 
SHOWN  FOR  COMPARISON 
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FIGURE  7.  DETAILED  COMPARISON  OF  TYPICAL  INTEKMITTENCY 
VELOCITY ,  AND  TEMPERATURE  PROFILES 


bursts  appear  at  the  edges  of  the  wake  (especially  it  high  x)  where  the 
sensible  axial  wake  motion  has  all  but  vanished.  In  Figure  7  for  example 
and  at  T)  ■  3.28,  u  ■  0.01  while  y  »  0.17.  We  shall  later  discuss  this 
phenomenon  further. 

In  order  to  proceed  with  the  computation  of  Y  and  J  which  involves  moments 
of  y,  the  distributions  of  Figure  6  were  curve-fitted  by  polynomials  of 
the  form 

V  -  £  AYn  (5) 

n  n 

with  the  aid  of  program  WEB-VI  written  for  the  Philco  2000  digital  computer 
(see  Appendix  A).  \dequate  fits  were  obtained  for  6th  to  8th  order  poly¬ 
nomials,  depending  on  the  XSTATION,  the  aim  being  to  minimize  the  computed 
standard  deviation  between  the  polynomial  and  the  experimental  curve. 

Table  I  shows  the  coefficients  A  for  each  XSTATION.  By  means  of  these 

n 

polynomials  and  the  aid  of  Equations  (2)  and  (3)  the  following  quantities 
were  computed: 

The  quantity  Y  was  computed  as  shown  in  Equation  (2)  and  is  pictured  on 

Figures  8  and  9.  The  former  shows  that  the  variation  Y  *  Y  (x)  follows 

closely  an  l/3-power  dependence,  which  is  natural  to  expect  because  the 

wake  itself  grew  as  x^.  The  Y  should  therefore  scale  as  L,  the  transverse 

—1/3 

wake  scale  which  grows  as  x  also  (Reference  9)  and  this  is  shown  on 
Figure  9.  The  average  of  the  points  lies  at  Y/L  *  2.11,  that  is  about 
two-th'rds*  of  the  wake  ’’radius"  out  from  the  axis  in  terms  of  the  non- 
dimei  clonal  radius  =  Y/L.  Figure  10  shows  that  Y  lies  very  close  to  the 
maximum  Y'  of  the  probability  distribution  dv/dY  implying  that  the  skewness 
of  this  distribution  is  negligible.  The  mean  standard  deviation  a  of  the 
front  corrugations  was  computed  from  Equation  (3).  Figure  8  shows  that  in 


♦Rather  than  three- fourths  as  Townsend  remarked  qualitatively 
(Reference  2,  p.  693). 
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FIGURE  8.  AXIAL  VARIATION  OF  FRONT  LOCATION  AND 
STANDARD  DEVIATION 
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FIGURE  9.  RELATIVE  MAGNITUDES  OF  THE  FRONT  LOCATION 
STANDARD  DEVIATION,  AND  TRANSVERSE  SCALE 
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FIGURE  10.  RELATIVE  MAGNITUDES  OF  THE  FRONT  LOCATION  AND 
THE  MAXIMUM  OF  THE  GAUSSIAN  DERIVATIVE 


dimensional  form  the  standard  deviation  also  follows  the  1/3  increase 
with  distance  past  about  x  =  40,  and  the  appropriate  scaling  with  L  is 
shown  in  Figure  9;  the  a  is,  on  the  average,  0.87L. 

Using  the  computed  value  of  c  it  is  possible  to  plot  the  distribution  of 

the  factor  y  about  the  front  position  Y,  as  shown  on  Figure  11.  With 

the  exception  of  the  initial  two  XSTATIONS,  all  data  points  fall  on  a 

well  defined  curve  which  is  in  fact  identical  to  that  encountered  for 
most  other  free  turbulent  flows  (Reference  1).  The  distribution  of 
intermittency  across  turbulent  interface  appears  to  universally  obey  the 
Gaussian  distribution  shown  in  the  latter  illustration.  The  agreement 
between  the  present  experiment  and  those  of  other  authors  will  be  shown 
in  a  later  illustration. 

5.2  MEASUREMENT  OF  ZEROS;  FRONT  MICROSCALE 

Figure  4  shows  the  circuit  used  to  read  the  frequency  of  zeros.  As 

expected,  very  few  events  are  noted  at  the  extremities  of  the  wake  or 

near  the  axis.  The  frequency  is  thus  zero  on  the  axis  and  outside  the 

wake  and  it  is  expected  to  peak  near  the  turbulent  front.  The  data 

are  shown  on  Figure  12  for  all  XSTATIONS  examined.  The  distribution  of 

the  zeros  N  about  the  mean  is  quite  symmetric  and  the  shift  of  the  Ymax 

(N  )  with  x  is  quite  clear.  The  peak  frequency,  originally  lying  at 
max 

about  20  kcps  gradually  shifts  to  about  30  kcps  far  downstream.  The 

and  Y  were  found  by  again  curve- fitting  the  data  of  Figure  12  with 

max 

polynomials  of  the  type 

cm"  (6) 

n 

n 

whose  coefficients  are  given  in  Table  II.  The  proper  polynomial  at  each 
XSTATION  was  again  found  by  minimizing  its  standard  deviation  from  the 
experimental  curve. 
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FIGURE  II.  DISTRIBUTION  OF  INTERMITTENCY  FACTOR  ABOUT 
FRONT  LOCATION 


(KCPS) 


COEFFICIENTS  FOR  CURVE-FITTING  POLYNOMIALS  FOR  ZERO  OCCURRENCES 
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-81.3  3.73  x  10  -9.01  x  10  1.16  x  10  -6.63 


Several  features  of  these  data  demand  further  examination.  It  is  natural 

to  suppose  that  the  gradual  spread  of  the  curves  of  Figure  12  with  x 

represent  the  physicaLly  enlarging  intermittent  zone  (i.e.,  the  standard 

deviation  c)  and  that  the  shift  of  the  position  of  N  corresponds  to 
_  r  max 

the  increasing  Y  along  the  wake,  as  shown  on  Figure  8.  In  studying  the 

latter  point,  the  ratio  Y/Y  was  plotted  as  shown  on  Figure  13  and  was 

found  to  agree  with  unity  to  within  4  percent  Furthermore,  the  frequency 

N  of  zero  occurrences  at  Y  (obtained  by  computing  N  from  Figure  12  at 
o  _  o 

the  position  Y)  is  identical  with  N  ,  as  shown  in  Figure  14.  The 
turbulent  front  position  corresponds,  therefore,  with  the  radial  distance 
from  the  axis  where  the  frequency  of  zero  occurrences  is  a  maximum. 

Finally,  the  symmetry  of  these  occurrences  was  checked  by  plotting  the 
normalized  zero  frequency  N/N  about  Y/Y  .  The  symmetry  of  these 
results,  shown  inFigure  15,  is  gratifying  and  the  resemblance  to  a  Gaussian 
distribution  is  quite  close. 

It  should  be  mentioned  that  the  measurement  of  the  zeros  in  the  manner 
described  is  neither  versatile  nor  unambiguous.  Much  depends  on  that 
portion  of  the  "intermi ttency  circuit"  which  rectifies  the  hot-wire  out¬ 
put.  If  the  rectifying  ( "mean- squaring")  process  is  too  slow  then  the  high- 
frequency  events  in  the  wake  (in  this  case  extending  to  some  hundreds 
of  kilocycles)  cannot  be  detected.  Of  the  process  is  too  fast  the  Schmidt 
trigger  fires  many  times  while  the  wire  is  within  the  turbulence,  not 
just  when  it  enters  or  leaves  it;  an  artificial  increase  of  N  is  thus 
observed.  By  contrast  the  latter  problem  would  not  affect  greatly  the 
measurement  of  intermittency  factor,  where  only  the  total  time  of  immersion 
in  turbulence  i.s  of  importance.  As  before,  the  storage  oscilloscope  was  of 
great  assistance  in  resolving  this  question. 

Because  both  the  intermittency  factor  y  3nd  the  zeros  N  are  distributed  in 
a  normal  (Gaussian)  fashion  as  shown  by  the  preceding  figures,  it  is  possible 
to  evaluate  the  autocorrelation  microscale  of  the  turbulent  front  in  the 
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FRONT  LOCATION  AND  THE 
ZERO  FREQUENCY 


FIGURE  14.  COINCIDENCE  OF  THE  ZERO  FREQUENCY  NQ  AT 
THE  FkONT  POSITION  AND  MAXIMUM  FREQUENCY 
N_  AT  EACH  X 
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manner  suggested  by  Corrsin  (Reference  1);  the  results  are  shown  on 
Figure  16,  normalized  alternately  with  the  local  transverse  scale  L  and 
the  local  wake  radius  Y.  It  will  be  noted  that  the  scaling  of  \  with 
these  two  parameters  (which,  of  course,  vary  as  x  ,  according  to 
previous  statements  and  the  results  of  Reference  9)  is  at  most  tenuous, 
although  some  assemblance  of  such  a  proportionality  appears  towards  the 
far  end  of  the  range  investigated  (x  -  90).  A  slight  but  still  inadequate 
improvement  is  noted  if  u(Y)  (i.e.,  the  mean  flow  velocity  at  the  front 
position)  is  used*  in  Equation  (4)  in  place  of  u#.  In  fact  it  would  seem 
from  Figure  17  that  v  is  probably  constant  along  the  wake,  as  illustrated 

r  _ 

by  comparing  \  with  the  (constant)  *Vc_,A. 

r  D 

The  apparent  lack  of  an  x-dependence  in  )  is  surprising  because  of  the 

F  -1/3 

strength  of  dimensional  arguments  relating  it  to  Ir^x  (all  other 
characteristic  lengths  such  as  A^,  A  ,  c,  V  etc.,  scale  with  L)  and  to 
a  lesser  extent  because  such  a  dependence  is  apparently  evident  in  bal¬ 
listic-range  experiments  (References  5,  6,  and  7).  One  possible  explana¬ 
tion  is  that  the  front  microscale  is  undergoing  a  relaxation  process  for 
x  <  90,  of  the  type  noted  by  Schapker  (Reference  5),  during  which  it 
remains  constant. 

In  this  connection  it  is  relevant  to  ask  whe<-Ui-r  the  scale  A  =  u  /N  itself 

F  ®  o 

does  not  have  a  physical  significance,  especially  because  of  the  repeatability 
observed  in  measuring  .  Since  Nq  is  the  "average  number  of  zeros"  A^  must 
have  the  meaning  of  an  "average  wavelength"  which  is,  of  course,  n//2  times 

larger  than  \  (cf  Equat ion  (4) ) .  Almost  exactly  the  same  numerical 

r 

difference  between  A  and  \  was  also  observed  by  Corrsin  and  Kistler  who 

r  c 

computed  from  measurements  of  the  average  extent  of  the  laminar  and 

*The  rigor  for  using  u#  in  Equation  (4)  rests  on  Corrsin  and  Kistler's 
conclusion  chat  "the  mean  velocity  everywhere  in  the  potential  part  of 
the  flow  must  be  constant  and  equal  to  that  at  infinity"  (Reference  1,  p.15). 
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turbulent  regions  at  Y  =  Y  (Reference  1,  p.  26  to  28).  This  disparity  was 

of  seemingly  minor  importance  in  the  latter  instance  because  no  additional 

phenomena  were  observed  to  give  A  a  physical  foundation  (in  the  same 

manner,  for  example,  as  connecting  the  microscale  to  the  Schmidt  trigger 

output).  In  the  present  case  additional  phenomena  featuring  a  wavelength 

A  were  observed;  these  will  be  described  farther  below. 

F 

6.  DISCUSSION 

6.1  EFFECT  OF  INTERMITTENCY  ON  FLUCTUATION  MAGNITUDES 

As  a  result  of  the  preceding  discussion  and  the  generally  accepted  view  of 

the  turbulent  front  structure,  the  output  of  a  stationary  sensor  in  the 
intermittent  region  can  be  thought  to  consist  of  the  bonafide  turbulent 
signal  superimposed  on  a  "clipped"  random  signal  due  to  the  in termi ttency , 
which  we  shall  henceforth  call  "pseudoturbulence".  Much  added  effort  is 
needed  to  study  each  of  these  two  phenomena  free  from  the  influence  of  the 
other.  In  the  meantime  we  shall,  in  this  Section,  examine  phenomena 

indicating  the  presence  of  pseudoturbulence  and  use  them  to  indicate  their 

effect  cn  the  fl’ctuation  magnitudes. 

2 

Asshown  in  Appendix  B  the  mean-square  (af)  of  a  variable  which  jumps  by 
H  across  the  front  is 

W)2  -  h\(1  -  Y)  (7) 

so  that  this  amount  is  added  to  the  overall  hot-wire  signal  due  to  the 
intermittency  (pseudoturbulence)  alone.  At  the  same  time  the  overall 
signal  also  should  be  reduced  by  the  familiar  factor  inversely  proportional 
toy  to  account  for  signal  "loss"  in  the  quiescent  zones  between  turbulence 
zones.  Psuedoturbulent  effects  for  which  the  wire  signal  needs  corrections 
consist  therefore  of  two  terms  of  opposing  sign,  as  first  brought  up  by 
Corrsin  and  Kistler  (eg.,  (57)  of  Reference  1). 
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We  now  ask  if  the  action  predicted  by  Equation  (7)  is  in  any  way  visible 
in  the  present  experimtnt.  Although  H  is  not  known  (and  most  likely 
depends  on  radial  distance)  the  factor  y(l  -  Y)  peaks  at  v  =  0.5,  i.e., 
at  the  front  position  Y.  Another  location  where  pseudoturbulent  effects 
are  visible  should  be  the  very  edge  of  the  wake  (T)  >  3)  where,  according 
to  Figure  7,  the  mean  velocity  and  temperature  have  acquired  the  stream 
levels  and  where  the  front  corrugations  are  too  narrow  to  contain  too 
many  of  the  small  eddies.  To  make  the  observations  at  these  two  loca¬ 
tions  easier,  we  recall  that  in  compressible  and  heated  flows  the  tempera¬ 
ture  and  density  jumps  (i.e.,  the  H  values)  across  the  front  are  certainly 
bigger  than  those  of  the  velocity  (Reference  1 3).  The  proper  observable  is 
therefore  the  temperature  -  in  fact,  it  is  the  ratio  (-T/T) / (v - 1 )M2 (lu/u) 
of  the  r.m.s.  temperature  to  r.m.s.  velocity  f luctuations*proper ly 
normalized  with  the  local  Mach  number  (Reference  IQ).  This  quantity, 
plotted  on  Figure  18  versus  radius  for  various  axial  positions,  should 
be  un.tv  in  a  homogeneous  turbulent  gas  without  heat- transfer .  In  the 
figure  both  the  maximum  of  v(l-y)  at  T)  V  2  and  the  precipitous  increase 
of  (AT/ T)  at  large  are  shown.  In  the  latter  case,  (Figure  18  (b))  we 
see  a  mechanism  driving  especially  large  temperature  fluctuations  near 
the  wake  edge  and  it  is  submitted  that  these  fluctuations  are  due  to 
p*eudoturbulence  as  anticipated  above. 

6.2  EFFECT  OF  INTERMITTENCY  ON  FLUCTUATION  SPECTRUM 

In  addition  to  the  effect  on  the  fluctuations  magnitude,  effects  on  their 
spectra  suspected  to  be  due  to  pseudoturbulence  have  been  observed. 

As  indicated  in  Reference  10,  the  longitudinal  spectra  of  the  turbulent 

fluctuations  of  wake  density  and  axial  velocity  have  exhibited  a  marked 

intensity  prominence  at  frequencies  ranging  from  about  30  to  50  kilocycles 

per  second.  It  was  concluded,  in  that  reference,  that  this  "peak"  in t  he 

spectra  was  due  to  the  front  corrugations  being  organized  into  a  weakly 

periodic  structure.  The  periodicity  was  found  to  be  most  intense  immediately 
*here  y  is  the  ratio  of  specific  heats. 
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FIGURE  18.  TEST  OF  THE  ’’STRONG  REYNOLDS  ANALOGY"  FOR  THE  RELAXING  ('TOP) 
AND  THE  SELF-PRESERVING  WAKE  (BOTTOM). 
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after  transition  to  turbulence,  decaying  o  that  it  was  barely  discernible 
at  X  *  90.  If  this  periodicity  was  due  to  the  pseudoturbu lence  it  is 
natural  to  expect  that  it  would  surface  clearly  if  the  pseudotui  oulerit 
spectrum  could  be  measured  separately  from  that  of  the  overall  signal. 

The  separation  of  the  two  types  of  spectra  is  not  an  easy  matter,  but  a 
fair  picture  of  the  situation  is  given  by  the  output  spectrum  of  the  Schmidt 
trigger  (i.e.,  of  the  lower  traces  in  the  photographs  of  Figure  5).  Figure 
19  shows  raw  data  of  the  trigger  output  spectra  compared  with  the  overall 
hot-wire  spectra  (i.e.,  the  top  traces  in  Figure  5).  Both  spectra  show 
clearly  the  peak  located  at  the  same  frequency. 

These  observations  suggest  that  the  front  is  organized  into  a  periodic 
structure,  and  can  be  further  discussed  with  the  aid  of  a  rudimentary 
model,  as  done  previously  by  Bradbury  (Reference  14)  and  as  shown  on 
Figure  20.  For  purposes  of  illustration  the  irregular  shape  of  the 
corrugated  front  is  replaced  by  a  periodic  structure  in  the  form  of 
adjoining  half-circles.  The  wire  signal  in  space  or  time  then  consists 
of  a  series  of  rectangular  waves  with  superpose j  turbulent  fluctuations; 
the  latter  is  very  much  in  accord  with  the  observations  (see  Figure  5) 
although  in  the  present  discussion  we  imagine  thac  the  rectangular  waves 
are  also  periodic.  It  is  natural  to  suppose,  as  shown  inFigure  20,  that 
the  spectrum  of  such  a  signal  is  somehow  decomposable  into  the  spectrum 
of  the  turbulence  itself  and  that  of  the  rectangular  waves  produced  by 
the  corrugations.  The  turbulence  spectrum  be  what  it  may,  the  pseudo- 
turbulent  spectrum  due  to  the  rectangular  waves,  discussed  in  Appendix  B, 
has  a  fundamental  at  the  characteristic  front  wavelength  and  weakening 
harmonics  at  integral  values  of  the  fundamental.  As  show.,  in  Figure  20, 
therefore,  the  spectrum  combining  the  turbulence  and  in termi t tency 
phenomena  should  show  a  peak  for  a  large  portion  of  the  radius  between 
the  axis  and  the  wake  boundary.  This  is,  in  fact,  observed  in  Figures  42 
through  63,  in  Reference  10.  Naturally,  the  observed  peak  due  to  the 
pseudoturbulence  is  "weak"  in  the  sense  that  its  Fourier  components  are 
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FIGURE  19  ,  QUALITATIVE  COMPARISON  OF  THE  SCHMIDT  TRIGGER 

OUTPUT  SPECTRA  WITH  THE  HOT-WIRE  OUTPUT  SPECTRA 
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FIGURE  19.  (CONTINUED) 


F09819  U 


FIGURE  20.  EFFECT  OF  INTERMITTENCY  ON  THE  SPECTRUM  OF  A  TURBULENT  WAKE 


in  reality  dispersed  cons iderab ly .  In  fact  the  Schmidt  trigger  output 
spectrum  of  Figure  19  looks  much  like  the  "Poisson"  spectrum  observed 
at  the  same  output  by  Corrsin  (Reference  1)  but  with  a  peak  superimposed. 

We  conclude  that  the  periodic  structure  is  indeed  tenuous,  perhaps  appearing 
and  disappearing  cyclically  as  suggested  by  Townsend  (Reference  2)  and 
Grant  (Reference  15). 

If  the  fluctuation  concentration  (the  "peak")  in  the  hot-wire  output 
spectrum  is  thus  identifiable  with  a  periodic  component  of  the  pseudo¬ 
turbulence,  it  is  worthwhile  to  test  Equation  (7)  for  that  component 
alone,  just  as  (7)  was  tested  in  Figure  18  for  the  entire  pseudoturbu lent 
signal.  Figure  21  shows  the  density  fluctuation  intensity  variation  with 
radius  for  a  narrow  (1  keps)  fluctuation  passband  about  a  frequency  at 
35  keps,  that  is  in  the  immediate  vicinity  of  the  preudoturbu lence  peak. 

The  maximum  of  this  curve  is  very  close  to  ~  *  2  where  the  maximum  also 
occurs  in  the  •(!->)  curve.  For  comparison,  the  mean-square  density 
fluctuation  integrated  over  all  frequencies  is  also  plotted;  its  own 
maximum  occurs,  at  this  case,  around  ~  *  1.5,  that  is  closer  to  the 
radial  location  of  the  maximum  shear  region,  as  already  discussed  in 
Reference  10.  It  appears  to  some  degree,  at  least,  that  the  behavior 
of  pseudoturbul ent  "periodicity"  is  described  by  the  rough  argument 
given  previously. 

6.3  WAVELENGTH  OF  THE  TURBULENT  FRONT 

The  periodicity  detected  in  the  wire  signal  and  its  clipped  (Schmidt 
trigger)  counterpart  has  a  wavelength  u(Y)/f^.  where  f_  is  the  peak 
frequency  in  the  spectra  of  Figure  19.  We  have  also  detected  another 
wavelength  u(Y)/Nq  where  N  is  the  frequency  of  zeros  at  the  front  loca¬ 
tion;  this  is,  as  explained  in  Section  5.2,  n// 2  times  the  front  micro- 
sale.  We  ask  what  the  connection  between  these  two  wavelengths  is,  and 
what  is  their  magnitude  related  to  other  characteristic  lengths  in  the 
wake . 
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FIGURE  21.  RADIAL  VARIATION  OF  THE  35* KILOCYCLE  FOURIER 
COMPONENT  OF  THE  DENSITY  SPECTRUM  (FOURTH 
FROM  TOP)  5HOWINC  TENDENCY  TO  FOLLOW  TMF 
PSEUDOTURBULENT  (THIRD  FROM  TOP)  RATHER  THAN 
THE  INTECRATED  TURBULENT  INTENSITY  ( BOTTOM) 


Figure  22  plots  these  two  wavelengths,  which  are  seen  to  be  about  equal 

farther  along  the  wake;  this  equality  could  in  anyway  be  predicted  earlier 

since  Nq  and  f^  were  found  to  be  about  equal.  Nc  explanation  exists  for 

the  differences  below  X  ■  50,  however.  A  value  A  =  u(Y)/N  =  u(Y)/f  = 

F  o  T 

commentary  on  A  will  be  made  further  below. 
F 


9+2  seems  to  prevail,  ^dded 


6.4  COMPARISON  WITH  OTHER  EXPERIMENTS 

It  would  be  interesting  to  compare,  whenever  possible,  the  present  results 
with  those  obtained  in  other  intermittency  measurements.  In  low-speed 
flows,  for  example,  Corrsin  and  Kistler  (Reference  1)  have  observed  that 
Y  is  distributed  (Y  -  Y)/~  in  a  manner  independent  of  the  type  of  floxv 
(wake,  jet  or  boundary  layer).  This  "universal"  curve  agrees  well  with 
the  present  results,  as  shown  on  Figure  23.  This  agreement  is,  of  course, 
fully  expected  since  y,  c  and  Y  are  in  essence  moments  of  each  other.  It 
is  hardly  surprising,  therefore,  that  measurements  subject  to  added 
natural  aberrations,  such  as  the  ballistic-range  measurements  of 
Levensteins  and  Krumins  (Peference  6)  show  equally  good  agreement. 

It  is  much  more  significant  to  relate  the  intermittency  behavior  to  the 
transport  of  momentum  across  the  turbulent  flow.  Thus  the  spreading  rates 
of  Y  and  ct  as  Corrsin  demonstrated,  follow  the  dependence  of  the  flow 
scale  L  on  distance  X  and  are  thus  dependent  on  the  type  of  flow  at  hand. 

—1/3 

It  has  already  been  seen  in  Reference  9  that  in  the  present  wake  L~X  ’ 

—  -1/3 

and  in  Figures  8  and  9  of  this  report  the  relationship  c~Y~L*-X  is  fully 

—  —i  j  2 

verified;  in  Townsend's  cylinder  wake  (Reference  12)  c~Y~L~X  as  expected. 


Table  III  shows  a  numerical  comparison  of  the  present  results  with  those 

available  from  low-speed  experiments;  the  significance  of  the  transverse 

scale  L  is  clearly  demonstrated  in  this  table.  It  should  be  noted  that 

the  position  Y  of  the  front  and  the  turbulent  macroscales  A  (the 

u 

longitudinal  velocity  scale)  and  A^  (the  longitudinal  temperature  scale) 
and  perhaps  the  front  "macroscale"  A^,  seem  to  be  universally  constant  for 
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FIGURE  22.  WAVELENGTH  OF  THE  TURBULENT  FRONT  NORMALIZED  WITH  THE  TRANSVERSE 
SCALE ,  AS  DEDUCED  FROM  WO  INDEPENDENT  MEASUREMENTS 


3 

co 
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FICURE  23.  COMPARISON  OF  THE  INTERMITTENCY  FACTOR  RADIAL 
DISTRIBUTION  WITH  OTHER  EXPERIMENTS 


TABLE  III.  COMPARISON  WITH  OTHER  EXPERIMENTS 


Present 

Work 

Townsend 

Grant 

Bradbury 

(Wake ) 

(Wake) 

(Wake) 

Ll£_Q _ 

Y/L 

2.1 

2.1 

2.1 

2.0 

a/L 

0.87 

0.80 

0.80 

0.44 

a/Y 

0.41 

0.38 

0.38 

0.22 

A  /  Y 
u 

0.5* 

0.5 

0.5 

0.55 

A  /L 

u 

l* 

1.06 

1.06 

0.92*** 

Vat 

0.5-1** 

— 

— 

0 . 58*** 

VL 

4.6 

— 

— 

— 

a/l 

F 

9_f2 

— 

— 

«»  w  • 

T](u  -  ,01)/L 

3.28 

— 

— 

•  •  • 

TKY  -  .01)/L 

4.5 

•  m  • 

•  •  m 

m  m  m 

★bad  scatter  (from  Reference  10) 

**0.5  on  axis,  I  near  edge  (from  Reference  10) 
***Corrsin  &  Uberoi  data  (Reference  19) 
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different  types  of  flow;  the  standard  deviation  differs  for  the  jet  and, 
as  Townsend  (Reference  2)  points  out,  for  the  boundary  layer  as  well. 
This  difference  has  been  attributed  by  Gartshore  (Reference  16)  to 
differences  in  the  effective  turbulence  Reynolds  number  R^,  among  such 
flows.  His  exa' r  expression,  put  into  the  present  nomenclature  reads 


a  2.97 
L  "  J  Rt 


(8) 


which  gives  a/L  =  0.83  for  the  R^  =  12.8  found  in  the  present  experiment 
(Reference  9).  N  )te  that  a  slight  difference  between  Equation  (8)  and 
Gartshore's  formula  (p.  97  of  Reference  16)  exists  because  the  transverse 
scale  L  is  about  227.  smaller  than  Gartshore's  Lq  which  is  the  flow 
"half  radius". 


The  possibility  of  a  periodic  front  has  been  well  established  by  Grant 
(Reference  15)  and  Keffer  (Reference  17);  Townsend  (Reference  2)  has 
made  calculations  of  possible  front  indentations  based  on  stability 
arguments  and  finds  that  such  waves  have  a  length  6.3L.  Recent  experi¬ 
ments  in  a  highly  heated  jet  (Reference  18)  have  disclosed  a  similar 
periodicity  in  its  interface  with  a  wavelength  of  about  8L.  We  can 
therefore  conclude  that  (a)  a  weak  periodicity  of  the  front  structure 
seems  to  exist,  (b)  the  wavelength  lies  around  9L  and  (c)  the  periodicity 
affects  the  spectra  sensed  by  the  hot-wire. 

As  previously  noted,  no  stretches  of  ambient  ( irrotat iona 1)  fluid  was 
found  on  the  wake  axis,  where  Y  *  1  always.  This  is  in  contrast  with 
the  results  and  Gibson  (Reference  3)  and  Baldwin  (Reference  4),  and  a 
tentative  explanation  of  the  difference  can  be  given  by  combining 
Gartshore's  results  (Reference  16)  with  a  brief  critique  made  by 
Behrens  (Reference  20)  of  low-speed  wake  experiments.  Behrens  observed 
that  the  turbulent  Reynolds  number  R^,  for  axi-symroctric  bodies  is  not 
constant  but  varies  with  bluntness,  from  about  1  for  very  blunt  (say, 
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disks  normal  to  the  flow)  to  about  10  for  very  slender  bodies;  at  the 
latter  extreme  R  is  about  equal  to  that  for  two-dimensional  wakes  (in 
the  present  wake  R  *  12.8  was  obtained,  very  close  to  the  value  of  12.7 
for  cylinders  at  low  speeds).  Further  along  the  scale  one  can  obtain 
R^,  as  high  as  30,  40  or  50  for  jets.  As  illustrated  on  Figure  24, 
Equation  (8)  predicts  large  differences  in  the  standard  deviation  Z  for 
flows  of  very  different  R  .  On  one  end  of  the  scale,  the  standard  devia¬ 
tion  for  jets  is  small  compared  to  the  flow  vidth.  For  the  present 
experiment  (R^,  =  12.8)  -  has  increased  to  the  point  where  the  axis  barely 
escapes  an  occasional  brush  with  the  external  flow.  For  blunt  axi- 
symmetric  bodies,  however,  the  "wrinkle  amplitude"  represented  by  c 
should  become  so  large  (it  is  z/ L  •  3  at  R^,  *  1;  in  the  present  work 
we  have  a/L  *  0.37)  that  the  front  corrugation  depth  reaches  and  exceeds 
the  wake  axis.  In  the  extreme,  therefore,  the  wake  would  not  resemble  a 
straight  cylinder  with  rough  surface  but  rather  a  cylinder  with  grossly 
ondulating  axis.  Gibson's  own  photograph  of  his  wake  (Reference  3)  seems 
to  confirm  this  very  description. 


A  comparison  can  finally  be  made  with  the  hypervelocity  observations  in 
ranges  as  summarized  inFigure  11  of  Reference  6,  if  we  put  i  *  2Y  in 
Figure  9, 


a 

6 


0.20 


The  agreement  is  poor,  and  is  definitely  due  to  the  manner  in  which  6  is 
defined  and  measured  in  the  range  experiments  (Reference  5,6)  and  especially 
the  possibility  of  underestimating  C  in  the  latter.  This  under-estimation 
must  be  in  turn  due  to  the  inability  to  detect  the  front  "valleys",  as 
mentioned  in  the  Introduction. 


The  comparison  of  the  microscale  is  difficult  again  because  the  range 
experiments  have  no  way  of  relating  with  L  or  “yC^A  where  the  latter 
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FIGURE  24.  RELATION  BETWEEN  STANDARD  DEVIATION  AND  TRANSVERSE  SCALE  FOR  FREE 
TURBULENT  FLOWS  AS  A  FUNCTION  OF  THE  TURBULENT  REYNOLDS  NUMBER 
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reters  Co  the  wake  alone.  However,  some  semblance  of  Che  dependence 

—  1/3 

^  ~  X  is  seen  in  Figure  13,  Reference  6. 

7.  CONCLUSIONS 

We  draw  che  following  conclusions  from  che  resulcs  presenCed  above: 

(1)  A  fully  Curbulenc  region  exists  abouc  Che  axis  of  che  wake; 

Chis  region  is  very  narrow  in  excenc  (compared  co  che  wake 
diamecer)  in  che  sel f- preserving  region.  No  occurrence  of 
irroCaCional  flov/  on  che  axis,  of  che  Cvpe  observed  by 
Baldwin  and  Gibson,  was  seen. 

(2)  The  wake  fronc  Y  is  locaced  ac  T]  *  2.1;  the  wake  "diamecer" 

2Y  thus  is  abouc  4.2L  where  L  is  the  transverse  scale.  The 
standard  deviacion  of  the  front  corrugations  C  *  0.87L. 

These  resulcs  are  in  very  good  agreement  with  low-speed 
wake  resulcs  but  in  poor  agreement  with  ballistic  range 
resu Its. 

(3)  The  distribution  of  y  about  the  front  position  (i.e.,  with 
(Y  -  Y)/c)  is  Gaussian  in  character  and  identical  with  all 
similar  variations  reported  so  far,  apparently  including 
ballistic  range  (hypervelocity)  experiments. 

(4)  The  intermittency  "radius"  of  the  wike  is  located  at 

TJ (y  ■  0.01)  ~  4.5,  contrasted  with  T|  (T  »  0.01)  «  3.65 
and  TJ ( u  *  0.01)  *  3.28,  the  "temperature"  and  "velocity" 
radius  respectively. 

(5)  The  longitudinal  autocorrelation  microscale  \  of  the 

F 

front  does  not  seem  to  scale  well  with  L  in  the  range 

of  X  studied,  although  at  large  X,  \  -  4L. 

r 

(6)  A  weak  periodicity  in  the  front  geometry  was  noted  with  a 
wavelength  Ap  •  9L  +  2L.  In  the  spectral  sense  this  periodicity 
affects  the  spectra  of  turbulence  within  the  wake  itself. 

The  characteristic  lengths  of  the  wake  arp  shown  on  Figure  25. 
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WAKE 


FIGURE  25.  SCHEMATIC  DIACRAM  OF  THE  WAKE  MODEL  COMPOSED 
FROM  THE  PRESENT  EXPERIMENTAL  RESULTS 
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APPENDIX  A 


THE  WEB-V’  COMPUTER  PROGRAM 


This  progr an  is  designed  to  process  intermi t toncy  data  in  wake- like 
turbulent  flows.  Its  primary  inputs  are  the  distribution  of  intermi ttency 
factor  V  (gamma)  and  zcro-occurrenco  frequency  N,  with  secondary  inputs 
from  the  WEB-I  and  WEB-V  programs. 

Referring  to  the  procedural  diagram  of  Figure  A- l,  a  number  of  discrete 

points  V  at  each  Y  are  inputed  at  each  XSTATION.  A  polynomial  y  (Y  ) 

Y  y 

is  produced  by  the  computer  to  curve-fit  this  distribution  (usually 
between  6th  and  8th  degree  polynomials).  The  curve-fitting  process  con¬ 
sists  of  searching  for  a  minimum  in  the  standard  deviation  between  the 
experimental  curve  and  each  polynomial;  the  polynomial  of  minimum  devia¬ 
tion  is  adopted  as  the  one  best  suiting  that  particular  XSTATION.  The 
average  front  position  Y  and  the  front  standard  deviation  cr  are  then 
computed  from  the  fitted  curves.  Various  non-dimensional  forms  of  Y 
and  a  are  finally  given  by  combining  with  other  wake  characteristics  as 
supplied  by  previous  diagnostics  of  the  axisymmetric  wake.  A  similar 
procedure  is  followed  with  the  table  of  numerical  inputs  (one  set  for 
each  XSTATION)  of  zero  frequency  N  versus  radial  position  Y^.  The 
necessary  curve- fitting  polynomials  are  again  at  most  of  the  8th  degree. 
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FIGURE  A- 1  WEB- VI  PROGRAM 


A- 2 


Q  k 


The  output  sheets  for  y  Include  the  following:  axial  positton  X  (XBAR), 

Y  (YBAR),  0  (SICMA),  ■  r/Y  (SIGMA1),  ^  -  c/L  (S1CMA2) ,  0 3  -  0/fC~A 
(SIGMA3)  ■  Y/L  (YBARi),  Y2  -  Y/ VcdA  (YBAR2)  and  Y',  which  is  the 
maximum  of  the  probability  curve  dr/dY;  also,  Y1  -  Y’/L  (Y1PR1ME) 

Y^  ■  Y’/  r CpA  (Y2PRIME),  and  the  fitted  curve  coefficients  and  standard 
deviation.  In  columnar  form  one  obtains  the  following  Y  (Y- GAMMA): 

Y  (GAMMA),  (Y-Y)Aj  (ZETA)  and  7)  (ETA). 


The  output  sheets  for  N  include  the  following  at  each  X:  Nq  (NO), 

Xfc  ■/ 2/nNo  ( LAMBDA- T) ,  Xj  -  X£  uffl  (LAMBDA-1),  Xj  -  >  ^  Vc^A  (LAMBDA-1), 
Xj  ■  Xj/D  where  D  is  the  cylinder  model  diameter  (LAMBDA-2)  X^  -  f.^/L 

/  *  i  itn  nt  1  \  \  _  /ii  v  s  «  i  _  •  -  .  .  i  _  _  -  _  _  ...... 


X./c  (LAMBDA- 5),  N  (N-M/\X), 
i  max 


(LAMBDA-3),  -  Xj/Y  (LAMBDA-4),  X& 

Ym,x  ’  YN  <W  <»*»«>•  N'l  ’  Wo  <N1>-  Y3  *  Ymax/Y  (Y3>  “nd  aUo 
the  coefficients  of  the  fitted  polynomial  and  its  standard  deviation. 

In  columnar  form  at  each  X  one  obtains  the  following  as  a  function  of  Y 

N>  Y4  ■  VWY4)'  N2  '  N/Nm.,x  <V  ’nd  *  <HA>  • 


N' 


Typical  results  are  shown  on  Figures  A-2  and  A-3. 
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FIGURE  A-3.  SAMPLE  OUTPUT  OF  THE  WEB-Vl  PROGRAM:  ZERO -FREQUENCY  RESULTS 


APPENDIX  B 


BEHAVIOR  OF  A  CLIPPED  RANDOM  VARIABLE 


B.l.  BASIC  STATISTICS 

Consider  a  stationary  random  function  f  (t)  such  as  shewn  on  Figure  B.l. 
This  "clipped"  variable  closely  resembles  the  time  history  of  the  output 
of  a  sensor  held  fixed  at  some  distance  from  the  axis  of  turbulence, 
except  that  the  tops  of  the  jumps  in  f(t)  are  now  considered  flat.  This 
va.’cform  then  resembles  even  closer  the  Schmidt  trigger  output  of  the 
interm  .ttency  circuit  employed  in  this  work.  Define  by  f^  and  f^  the 
values  of  f  at  the  two  extremes  (related,  say,  to  the  sensor  being  out¬ 
side  and  inside  the  turbulence  respectively)  separated  by  a  constant  H. 
Similarly  ut^  and  refer  to  time  intervals  outside  and  inside  the 
turbulence,  and  intermittency  y  over  a  long  Interval  T  is  defined  by 


Y 


T  s  DAtj  +  EAt2 


and 


IAt2  T  - 
T  "  T 


l  -  V 


B-I 


The  moan  value  of  f(c)  now  is 

*  T  -Q  fdt  "  T  “flutl  *  T  Lf2“C2 

-  (1-V)fl  *  Vff^H)  ■  fL  4-  HV  (B-l) 

which  reduces  to  the  obvious  results 

f <V-0)  -  ti 

f(y-l)  -  ti  4-  H 

The  mean-square  of  the  variance  (fluctuation)  of  £(t)  is 

.  . .  j 

<Af)2  «  u-7]2  .if  <f-7)2dc 

1  o 

with  f  equal  to  f^  t  H  or  f^  as  appropriate  and  f  given  by  (B-l). 
The  result  is 

(if)2  -  YH2(L-\)  (B-2) 

Here  we  distinguish  two  important  limits.  First,  let 


which  gives 


which  is  infinite  at  y  ■  0.  Physically,  this  is  to  be  expected  at  the 
extreme  edge  of  the  turbulence  zone  for  situations  where  the  external 
(ambient)  value  of  the  observable  f  is  zero.  Such  a  case  should 
theoretically  occur,  for  example,  when  f  is  a  contaminant  density  and 
the  environment  is  completely  devoid  of  this  contaminant. 


B-3 


The  other  extreme  occurs  when  »  H. 


Then,  f  -  f^  and 


Vd-Y)  7-  I 
1 


This  value  Is  representative,  for  example,  In  the  present  wake  experiment 
where  the  Jumps  In  all  observables  (velocity,  temperature,  etc.)  are 
rather  weak  across  the  interlace. 


B.2.  PERIODIC  VARIABLE 

Equations  (D-l)  and  (B-2)  are  true  regardless  how  the  rectangular  waves 
of  f(t)  are  arranged.  Consider,  now,  that  the  pattern  repeats  It  '.elf  In 
time  at  regular  Intervals  2^/1  ot  circular  frequency  -A;  the  Intermit tency 
factor  Y  retains  Its  previous  definition.  Wc  are  interested  in  the 
spectrum  of  this  periodic  clipped  function.  Toe  Fourier  series  describing 
the  function  is  given  (Reference  21)  by 


f(t)  •  yh  +  ^ 


2H  - 

L. 

n 

n* 


sin  nr". 


tcostutcos 


n~v 

2 


-  sin  nit  sin 


rrv -i 
2  J 


■  YK  +  E  (A  cos  nit  +  B  sin  mt) 

^  n  n 

where  H  is  the  constant  height  of  the  rectangular  waves.  The  spectrum  of 
f(t)  can  now  be  constructed  since  the  mean  square  intensity  of  the  Fourier 
components  is  given  by 

2H^  sin^  nry  _ 

— - - 2 a  mean- square  of  component  n 

it  n 

which,  of  course,  becomes  zero  for  Y  ■  0,1.  In  this  spectrum  the  most 
intense  component  is  the  fundamental  (n  •  l)  with  gradually  weakening 
higher  harmonics  at  integral  values  of  the  fundamental  frequency. 


